Progress in the mathematical theory of 
quantum disordered systems 

Walter F. Wreszinski 
Instituto de Ffsica 
Universidade de Sao Paulo 
Caixa Postal 66318 
05314-970 Sao Paulo, SP , Brazil 

November 1, 2012 

Abstract 

We review recent progress in the mathematical theory of quantum 
disordered systems: the Anderson transition , including some joint 
work with Domingos Marchetti, the (quantum and classical) Edwards- 
Anderson (EA) spin-glass model and return to equilibrium for a class 
of spin-glass models, which includes the EA model initially in a very 
large transverse magnetic field. 

1 - Introduction 

In recent years there has been a significant progress in the mathematical 
theory of (quantum) disordered systems. Our purpose here is to present the 
main ideas (without proofs), with a clear discussion of the ideas, concepts 
and methods. Some new remarks and results pertaining to sections 3 and 4 
are also included. 

We shall be interested in properties of disordered systems at low tem- 
peratures, i.e., near the ground state (absolute temperature T = 0); in the 
critical region of spin glasses, there exist the spectacular recent rigorous re- 
sults on the mean field theory (see the review by F. Guerra |Gue06j ). and 
various results for other disordered systems, including the well understood 
high temperature phase in spin glasses, are discussed in the comprehensive 
book by A. Bovier [Bov06j. 



Our restriction to very low temperatures implies, of course, that we shall 
be dealing exclusively with quantum systems. In particular, the Ising model, 
when it appears, must be regarded as the anisotropic limit of quantum (spin) 
systems: the fact that this is not only mathematically so is emphatically 
demonstrated by the well-known fact that the critical exponents of the Ising 
model in three dimensions (see, e.g., [ZJ79] ) are surprisingly close to those 
measured in real magnetic systems, precisely because most of the latter are 
highly anisotropic. 

Three important issues appearing in the above-mentioned context are lo- 
calization, arising in connection with the Anderson transition, which we 
discuss in section 2, reporting on several works, as well as on some joint work 
with Domingos Marchetti ( [MW12a] . |MW12b] ) . frustration, which appears 
in short-range spin glasses such as the Edwards-Anderson (EA) model, dis- 
cussed in section 3, and the (thermodynamic version of) instability of the 
first and second kind, which relates to the return to equilibrium for special 
initial states and probability distributions in a class of models (including the 
EA model) in section 4. 

2a - The Anderson transition - Generalities and review of some 
results 

The breakdown of translation invariance (or, more generally, Galilean 
invariance in many-body systems with short-range forces) leads to the exis- 
tence of crystals, and the corresponding Goldstone excitations are phonons 
|Swi67j . On the other hand, the discrete translation invariance subgroup of 
the crystal is also frequently broken by e.g. impurities. This fact brings 
about a number of important new conceptual issues, and is usually modelled 
by the introduction of a random local potential in a tight-binding model (for 
the latter see |Fey65| ). The associated physical picture consists of lightly 
or heavily doped semiconductors (e.g. Si doped with a neighboring element 
which contributes excess electrons, for instance P). This is the Anderson 
model ([And58], |And72j ). described by the Hamiltonian 



A + V L 



(1) 



on l 2 (Z d ) where A is the (centered discrete Laplacian) 




n':\n— n'\=l 
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plus a perturbation by a random potential 



(2) 



(3) 

where {V^} neZ d ^ s a faniily of independent, identically distributed random 
variables (i.i.d.r.v.) on a probability space (Q,B,fi), with a common dis- 
tribution = fj,({u> : < x}); is assumed to depend linearly on a 
quantity v > 0, which is the disorder parameter, also called coupling con- 
stant. The spectrum of H w is, by the ergodic theorem |KS81] . almost surely 
a nonrandom set a(H w ) = [—2d, 2d] + supp dF . Anderson conjectured that 
there exists a critical coupling constant < v c < oo such that for v > v c 
the spectral measure of H is pure point (p.p) for /^-almost every u, while, 
for v < v c the spectral measure of H w contains two components, separated 
by so called mobility edge E . if E G [E~~,E + ] the spectrum of H 1 ^ is 
pure absolutely continuous (a.c); in the complementary set a(if a; )\[i?~, E + ], 
H u has pure point spectrum - leading to the new important phenomenon of 
localization, first proved for d = 1 in |GMP77j . and for d > 3 for the first 
time in by |JF85j , based on previous seminal work by [FS83] , for fixed energy 
and large disorder, or for fixed disorder and large energy (the latter meaning 
near the edges of the band, i.e, the boundary of a(H u )); see also [AM93] for 
a much simpler proof by an entirely different method, and references given 
there. 

The references given above are, of course, only symbolic for the extensive 
literature on the Anderson model, but, fortunately, an excellent review on 
the subject, up to about 2007, is available |Jit07] . to which we refer for a 
readable and enlightening discussion, as well as references. Our review on the 
Anderson model will be mostly restricted to later developments, related, in 
particular, to the elusive existence of absolutely continuous (a.c.) spectrum 
and the Anderson transition. Even with this restriction, our exposition will 
not be exhaustive, but we shall try to explain some of the key issues, with 
the aim of stimulating further research in this exciting area. 

The first model for which a sharp metal-insulator transition was found is 
the one described by the almost Mathieu operator H^xe acting on Z 2 (Z) and 
given by 

(# WiA ,0*)(n) = *(n + 1) + *(n - 1) + /(wn + 9)V(n) (0) 
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with f{6) = Acos(27r#). This model was extensively studied in the physics 
and mathematics literature since the seventies, and displays, for almost every 
to G R, 6 G R, a transition from a set of delocalized states ("metal"), 
characterized by a purely a.c. spectrum for A < 2, to a set of localized 
states ("insulator"), characterized by a purely pure point (p.p.) spectrum 
with exponentially decreasing eigenfunctions (see jJitOTj .3.4) for a rather 
complete review and references). The one-dimensional nature of the model, 
as well as its physical range of applications (two-dimensional electron subject 
to a perpendicular magnetic field, integer quantum Hall effect) distinguish it 
from the Anderson model (see, again, |Jit07| . pg. 632 and references). 

A first direction in the search of a.c. spectrum led to the study of models 
with randomly decaying potentials, where V£ in (3) is replaced by a n V^, 
with dfi y as | ti | y oo . For dimension d — 1, many sharp results were 
obtained, for d > 2, see |Krib] . |dMS99j . and the lecture note |Bou03] . For 
the special "random surface case", where V£ = outside a " hypersurface" 
(of some thickness) in Z d , see [JLOlj . In a different, but related, approach 
[KS09J studied the eigenvalue statistics of a class of matrices (CMV matri- 
ces, see |Dam07j ) with random decaying coefficients, and [Sch09], random 
band matrices, with entries vanishing outside a band of finite width around 
the diagonal, with a view to determining the crossover between an "insulat- 
ing" regime with localized eigenfunctions and weak eigenvalue correlations 
and a "metallic" regime with extended eigenfunctions and strong eigenvalue 
repulsion. 

Much more recently, Krishna jKrid] presented a class of Anderson type 
operators with independent, non- stationary (non-decaying) random poten- 
tials, with pure a.c. spectrum in the middle of the band for small disorder, 
and |Kricj a class of continuous, non-decaying and non-sparse potentials ex- 
hibiting a spectral transition. These papers also open interesting possibilities 
of applications of well-known methods, powerful in other contexts (positive 
commutators, virial theorem, Mourre theory) in the context of the Anderson 
transition. 

In the nest section 2b, we confine our attention to two special models of 
the Anderson transition. One is the Anderson model on the Bethe lattice, 
first sudied by A. Klein in a seminal paper [Kle98j (see also |Jit07j for a 
discussion and further references). This work was further considerably ex- 
tended by M. Aizenman and S. Warzel f jAWllaj . |AW06] . jAWllbj . [AWb] . 
| AWaj ) . The other is a sparse model in d > 2 dimensions, studied by D. 
Marchetti and the present author [M W12a] , |MW 12bj. Our special emphasis 
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on these two models is justified by the fact that they are amenable to a quite 
detailed analysis ; it does not, of course, imply any comparative judgment 
og value regarding the other approaches mentioned in section 2a. Although, 
as we shall see, both models have serious shortcomings when compared with 
the original proposal (l)-(3), both yield a pallette of stimulating insights. 
We treat the sparse model first, because it leads to the introduction of some 
notions which are used elsewhere in the sequel. 
2bl-Sparse models 

Another collection of models introduced to study the Anderson transition 
have V£ in (3) of the form 



with elementary potential (" bump") (p u : 7L d — > R satisfying a uniform 
integrability condition 



\^(z)\ < 



Co 



1 + \z\ d+£ 



for some e > and < C < oo and 

lim # { i:M<fl} =0 

R-^oo R d 

(4) 

Due to condition of zero concentration, potentials such as V are called sparse 
and have been intensively studied in recent years since the seminal work by 
Pearson |Pea78j in dimension d = 1 by Krishna [KriaJ, Krutikov [Kru04j and 
Molchanov [Mol98j, [Mol99] in the multidimensional case. For d > 2 the 
interaction between bumps is weak while for d — 1 the phase of the wave 
after propagation between distant bumps becomes "stochastic". We shall 
report on a version of these models recently studied by D. Marchetti and the 
present author |MW12a] . |MW12b] . for which a rather detailed information 
is available. We consider (infinite) Jacobi matrices 

(J u) n = u n+1 + u n _i , 

with the perturbation potential 
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with u = (u n ) n > G / 2 (Z+) ,and 




v if n = a'j G .4 , 
if otherwise , 



for v G (0, 1). .4 = {cij }j>i denotes a random set of natural numbers a" = 
dj + Wj- , where a, satisfy the "sparseness condition": 



with ai+1 = (3 > 2 and Uj, j > 1, are independent r.v. on a probability space 
(0,5, /i) uniformly distributed on a set Aj = {—j,. ■ ■ ,j}- These variables 
introduce an uncertainty in the positions of those points for which v n ^ 0: 
such models are called Poisson models (see |Jit07] and references given 
there). Note that the support of the Uj only grows linearly with the suffix j. 
We write J w = J + V" '■> an d denote by the operator associated to J w on 
the Hilbert space % of the square integrable sequences u = (u n ) n> _ 1 which 
satisfy a 0-b.c. at —1: 



The essential spectrum of equals [—2,2]: it will be represented as 
A = 2 cos a with a G [0, 7r). Zlatos |Zla04j proved that this model exhibits a 
sharp transition from s.c. to p.p. spectrum. This was shown independently 
in f |dCMWll| . [MW12bp : 

Theorem 1 Let J% be as above. Let 



with v G (0, 1) and (3 G N, f3 > 2 . Then, for almost all u with respect to the 
uniform product measure on A = {— j, . . . , j}, 

(a) there exists a set Ai of Lebesgue measure zero such that the spectrum 
restricted to the set I\Ai is purely singular continuous, 

(b) the spectrum of Jp^ is dense pure point when restricted to 

I c = ([— 2, 2]\2 cos 7rQ) \I for almost every G (0,7r), where characterizes 
the boundary condition. Thus it is purely p.p. in this interval. 
Remark 1 



dj — Oj_i = f3 3 , 



j = 2,3,... 



m_i cos — Mo sin = 



1= AG [-2,2]\2cosvrQ : — (B - 1)(4 - A 2 ) > 1 

11* 
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Figure 1: Figure 1: In the above figure, we see the s.c. (light gray) and p.p. 
(dark grey) spectra separated by the "mobility edges" X^ = ±2y(l — zr), 
with = 1.3038---. 

The basic emphasis on one-dimensional (sparse) models has an important 
technical reason: the profound approach of Gilbert and Pearson connecting 
the space asymptotics of eigenfunctions of the restrictions of a large class of 
Sturm-Liouville operators to finite intervals to the spectral theory of these 
same operators in infinite space through the concept of subordinacy [GP87J 
is only available in one dimension. Using the important transfer matrix 
version of this theory due to Last and Simon |LS99j . the surprise is that, 
even in a regime of strong sparsity, a spectral transition from s.c. to p.p. 
spectrum (first shown by Zlatos |Zla04] ) may be proved: this is theorem 1. 
Our proof (theorem 4.7 of |MW12b] ) differs from that in |Zla04j by the use 
of the (optimal) metrical version of Weyl's theorem on uniform distribution 
|KN74j (in this case, of the so-called Priifer angles, see [KLS98J) due to 
Davenport, Erdos and LeVeque [DEL63J. The latter permits an explicit 
characterization of the exceptional set in which the spectrum may not be 
of p.p. type (this exceptional set arises in connection with the concept of 
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essential support of a measure |GP87j ). which reveals it to be also of p.p. 
type. As a result, the spectrum at the edges of the band is proved to be 
purely p.p.. 

The robustness of the transition depicted in theorem 1 follows from |dRJLS 96j , 
because the Hausdorff dimension of the s.c. spectrum may be seen to be 
nonzero f jZlaOi] . jdCMWlOj l 

Finally, the difference Laplacean (2) leads to tunneling and bands [Fey65|, 
with a pure a.c. spectrum. Starting from the difference Laplacean, the 
figure shows that any strength of the disorder v, however slight, yields p.p. 
spectrum in a nonempty range of energies. In addition, this occurs in a 
regime of strong sparsity. This is also true for the forthcoming models in 
d > 2 dimensions, and is thus a concrete manifestation of the tunneling 
instability, first studied by Jona-Lasinio, Martinelli and Scoppola |LMS81j . 
and called by Simon a "flea on the elephant" [Sim85j, a phenomenon which 
has also dynamical counterparts [GSOlJ, [WC98J, and is at the root of the 
existence of the chiral superselection rules induced by the environment which 
account for the shape of molecules such as ammonia, see the review by A. S. 
Wightman [ Wig95| . 

Unfortunately the s.c. spectrum does not posess either the dynamic or the 
perturbation theoretic properties (see |SW86j . |How80j for the latter) which 
are commonly associated with the physical picture of delocalized states. Re- 
garding the dynamical properties, for instance, the sojourn time of a par- 
ticle with initial state \I> in a compact region S is defined as 

/oo 
\\P s exp(-itH)V\\ 2 dt (5) 
-oo 

where Ps denotes the projector on S. If \& G % sc , it follows from a theorem 
by K.B. Sinha [Sin77j that a S exists with \S\ < oo and J(S] = oo. See 
also |MW12aj . |MW12bj for further discussion. 

We try therefore to attain higher dimensions. 

Multidimensional version 

Consider the Kronecker sum as an operator on 1-1®%: 

4 2) ■= Jf ®i + ei® jf 

(6) 

where to 1 = (u;j) . >± and co 2 = (u;J) >1 are two independent sequences of 
independent random variables defined in (Q, B, u), as before (we omit u 1 and 
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oj 2 in the l.h.s above for brevity). Above, the parameter 9 e [0, 1] is included 
to avoid resonances . We ask for properties of (e.g. the spectral type) 
which hold for typical configurations, i.e., a.e. (a; 1 , u 2 , 6) with respect to 
v x v x I where I is the Lebesgue measure in [0, 1]. J e is a special two- 
dimensional analog of J^; if the latter was replaced by — A + V on L 2 (IR, dx) 
where A = d 2 /dx 2 is the second derivative operator, and V a multiplicative 
operator Vip(x) = V(x)ip(x) (potential), the sum above would correspond 
to (—d 2 /dx\ + Vi) + {—d 2 jdx\ + V 2 ) on L 2 (R 2 ,dxidx 2 ), i.e., the "separable 
case" in two dimensions. Accordingly, we shall also refer to jj} n \ n = 2, 3, . . ., 
as the separable case in n dimensions. 

Our approach was to look at the quantity 



(*,e-*^*) =f 1 {t)f{9t) 
where 

/*(«) = £(*)+&(*). i = l,2 

with&(a) = / e- isX d^(X) 



Above $, * en®n, 



with ^ e "H sc , pj, Xi € "Hpp and f+p denotes the direct sum of two vectors 

tp,peU. 

We shall use the following folklore proposition (see, e.g., theorem 3.2 of 
|MW12bj or [3m77j): 
Proposition 1 

Let p, be a measure on the space M(R) of all finite regular Borel measures 
on K. If the Fourier-Stieltjes transform of p, 

K 9 t 1 — ► /2(t) = / e- itx dfi(X) 
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belongs to L 2 (IR, dt), then /i is absolutely continuous with respect to Lebesgue 
measure. 

The time-like decay of the Fourier-Stieltjes (F.S.) transform of the spec- 
tral measure is dictated by the Hausdorff dimension of the spectral measure: 
Let 5* be a subset of R, a G [0, 1], and 5 > 0. Define 

oo 

\S c U^Cj with \Cj \ < 5} 
where \C\ denotes the Lebesgue measure (length) of C, and 

h a {S) = limh%(S) = suph%(S) 

We call h a a - dimensional Hausdorff measure on R. h l agrees with Lebesgue 
measure, and h° is the counting measure, so that {/i a |0 < a < 1} is a family 
which interpolates continuously between the counting measure and Lebesgue 
measure. 

Definition 1 A Borel measure on R is uniformly alpha-Hoelder con- 
tinuous (UaH) if there exists C such that for every interval / with |/| < 1, 

H{I) <C\I\ a 

By an important theorem of Last [Las96], UaH measures may be ob- 
tained by a process of closure. In our one-dimensional model the Hausdorff 
dimension varies locally in the s.c. spectrum, and the local Hausdorff dimen- 
sion (suitably defined) may be determined explicitly. There exists a dense set 
in the s.c. subspace such that the spectral corresponding spectral measure 
is UaH. We shall for simplicity assume that the Hausdorff dimension has a 
constant value which will be denoted by a. 

The basic theorem, due to Strichartz [Str90j (with a very slick alternative 
proof by Last |Las96j . therefore we call it Strichartz-Last theorem), relates 
decay in the Cesaro sense to the Hausdorff dimension: 

Strichartz-Last theorem 

Let fi be a finite UaH measure and, for each / e £ 2 (R, d[i), denote 
ffi(t) = \ exp(-ixt)f(x)dfj,(x) 
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Then there exists C depending only on fi such that for all / G L 2 (R, dfi) and 
T > 0, 

(\M 2 ) T < C\\f\\ 2 T-« 

where ||/|| denotes the L 2 - norm of /, and (g)x = 9 ^ dt . 
A measure is called a Rajchman measure iff 

lim ji(t) = 

\t\— toe 

It does not follow from the decay in the Strichartz-Last theorem that the cor- 
responding measure is Rajchman . Let E denote the usual Cantor "middle- 
thirds" set in [— 7i, n]] die F.S. transform of the corresponding measure F 
is: 

oo 

T(u) = JJcos[2/3mtt3^ +1 ] 

i=i 

The corresponding Hausdorff dimension is well-known to be a = j jpg^j und 
the measure is UaH, but, from the above formula for T it follows that 

r(n) = r(3n) for all n G Z 

and hence the Cantor measure is not Rajchman. 
The main result of |MW12a] was: 
Theorem 2 Let 



v 2 < a (v 7 /?- l) < v 2 c 



with a < 4, where v c = 1\flj5 — 1). Then, for almost every (a; 1 , a; 2 , 9) with 
respect to v x v x I, [a.] there exist A with A + = — A~ and 

< A + < A+ 

such that 

(A-(l + fl),A + (l + #)) Ca ac (jf) 

[b.] 

[-2(1 + 9), A-(l + 9)) U (A+(l + 9), 2(1 + 0)] c a pp (jf 
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[a] 



<x sc (jf)n(A-(i + #),A+(i + #)) 

may, or may not, be an empty set. 

The basic idea of the proof is that the Strichartz-Last theorem suggests a 
pointwise decay of the F.S. transform of the spectral measure of type t~ a ^ 2 for 
large t. For the Kronecker sum (with 6 = 1) the F.S. transform is the product 
of the corresponding transforms for the one- dimensional system, which we 
expect decays as t~ a . In order to use the proposition, we restrict the spectral 
measure to an interval in the s.c. spectrum such that 2a > 1 which is in 
principle non-empty (and may be proven so), and denote the F.S. transform 
of the thus restricted measure by the same symbol as before. It turns, how- 
ever, out that this heuristics is not correct mathematically: as shown above, 
Cesaro decay does not in general imply pointwise decay or, in other words, 
that the spectral measure is Rajchman. This is a much harder problem - (see 
|MW12b] . chapter 4 for the full treatment of the pointwise decay of a model 
with superexponential sparsity and pure s.c. spectrum, and the recent very 
hard analysis to prove local decay in nonrelativistic QED |TJFS12j ). More- 
over, the result for general 9 may not be true due to resonance between 
Cantor sets, a subtle phenomenon which has been so far analysed only in 
the self-similar case by Hu and S.G. Taylor [HT94j. Self-similarity occurs, 
however, seldom: indeed, the s.c. spectrum of sparse Jacobi matrices is 
not self-similar by a theorem of Combes and Mantica |CM01j . We proved, 
however, that the main idea is indeed correct by generalizing a method due 
to Kahane and Salem ( |KS94j . |KS58j ) in specific (self-similar) cases to the 
problem at hand. 

We shall come back to this model in section 2d. 

2b2 - Anderson model on the Bethe lattice 

Consider a particle moving - for simplicity - on the Bethe lattice, a regular 
rooted tree graph, with vertex set Q, in which each vertex, aside from the 
root O has K + 1 neighbors with K > 2, and an edge set £, described by a 
Schrodinger operator of the form 

(H x (u)*)(x) = - *(v) + W u (x)*(x) (7) 

(x,y)e£ 

V u is a real valued random potential, with u again representing the ran- 
domness. We assume, again for simplicity, that the {V u (x)} xe g comprise 
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i.i.d. random variables, whose probability distribution is of the form p(V)dV 
with p the uniform distribution, constant over the interval [—1,1] and zero 
elsewhere (the uniform distribution). The regularity conditions assumed in 
|AW1 lb] . |AWb] . |AWa] are much more general, but the uniform distribu- 
tion is always a possible special case, and is the simplest example to be kept 
in mind. H\(u) is a bounded self-adjoint operator on l 2 (Q), with spectrum 
which is, for almost all realizations of u, the nonrandom set 

a(H x ) = [-\E X \, \E X \] with E x = ~{2^{K) + A) (8) 

The fact that the spectral measures associated to the different vectors 5 X G 
l 2 (calG), where 5 X denotes the Kronecker delta at x, are almost surely all 
equivalent, i.e., a.s. cyclicity of the spectrum, is an important result of plOOj. 

For a tree with N generations, there are N$ = (K + VjK 1 ^' 1 vertices at 
the surface (generation JV = 0) and a total N T = (K N+1 + K N ' 2 )/(K - 1) 
vertices. Since the ratio N$/Nt is nonzero for N — > oo, the behavior of 
statistical systems on a Bethe lattice is quite distinct from those on a regular 
Bravais lattice. 

Models such as (7) have been considered since the early days of Anderson 
localization, see [And58j , [ RACT73] . |ACT74j . The first mathematical proof 
of a transition was the seminal result of A. Klein |Kle98] : 

Theorem 3 (Theorem 1.1 of |Kle98j ) For any E, < E < 2xf{K), there 
exists X(E) > such that for any A with |A| < X(E) the spectrum of H\ in 
[-E, E] is purely a.c. with probability one. 

Aizenman |Aiz94j proved localization in the Bethe lattice for energies 
beyond K + 1, at weak disorder: this quantity appears as the edge of the 
I 1 - spectrum of the free Laplacean on the lattice. What happens between 
2^{K) and K + 1 for small A? (This is problem 4 in Jitomirskaya's list of 
open problems). As noted in |Jit07] . bottom of page 626, in analogy with the 
random surface case or the random decaying potentials mentioned in section 
2a, one might conjecture that l l - phenomena prevail up to the I 2 - threshold, 
i.e., that there is localization up to near the edge of the I 2 spectrum - but see 
the following theorem 4b, due to |AWllbj . |AWb] . a surprising result which 
contradicted original expectations, and whose ideas we now briefly sketch. 

As explained in |AWllaj . |AWllbj . and stated heuristically, states which 
may locally appear to be localized have arbitrarily close energy gaps AE 
with other states to which the tunneing amplitudes decay exponentially in 
the distance R, as exp(— L\(E)R), where L\(E) is the Lyapunov exponent 
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defined as follows. Let 

G x (x,y;(,uj) = (5 x ,(H x (uj) -()-%)) 
where £ G C + = ( G C s.t. 3^ > denote the Green's function. Then 

L X (E) = -\imAv(\log\G x (0,0;E + ir],uj)\) (9 - a) 

T)\0 

where Av stands for the expectation over the u> probability distribution. 
Indeed, it is shown in |AWb] that on tree graphs L x provide the typical 
decay rate of the Green's function: 

\imG x (0,x;E + iri,oj) ~ C exp(-L x {E)\x\) (9-6) 

a.s. in u, where the limit exists for a.e. energy E and \x\ = d(x,0) denotes 
the graph distance of the vertex x to some fixed but arbitrary vertex O G Q . 
Mixing between two levels would occur if AE <C exp(— L X (e)R), which is 
a very stringent condition, not satisfied in finite dimension, but, since, as 
previously remarked, the volume grows exponentially, as K R = exp(i?log K), 
extended states may form, and a sufficient condition may be shown to be 
L x (E)>\ogK: 

Theorem 4a (Theorem 1 of |AWllb] , see also |AWb] Under the previous 
assumptions, at a.e. E G cr{H x ) at which also 

L X {E)> log K (10 -a) 

one has, with probability one, 

SG x (x,x;E + i0,u) > (10-6) 

for all x G Q. 

Above, E+iO denotes lim^o of the corresponding quantity with argument 
E + ir]. It is well known that 

ir- l %G x {x, x; E + iO, to) (10 - c) 

is the density of the a.e. component of the spectral measure associated with 
the state \x >= 5 X (see, e.g., |Jak06] . Theorem 3.15), and thus theorem 4a 
implies that in any energy interval within a(H x (u)) for which (10a) holds, 
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the random operator a.s. has some a.c. spectrum on a subset of energies of 
positive Lebesgue measure. We now have: 
Theorem 4b (Theorem 2 of |AWllb] ) For 

A < A K = (yf{K) - l) 2 /2 
condition (10-a) holds for all energies satisfying 

\E\>\E x \-5{\) (10 -e) 

at some 5(X) > 0. 

Clearly (10-e) includes the spectral edges [E x , E x +S(X))u(\E x \-S(X), \E X \]. 
Thus, for A small enough (as in theorem 4b), near the spectrum's boundary 
the random operator has a.s. only purely a.c. spectrum, and thus there 
is no mobility edge beyond which the spectrum is pure point and local- 
ization sets in - a surprising result, which contradicts original expectations 
|RACT73][XCT74] . On the basis of the joint continuity of L X (E) in (A, E), 
it has also been conjectured |AWllb] that for small enough A - for bounded 
potentials and sufficiently regular probability distributions - the spectrum is 
purely a.c. 

An improved version of theorem 4a includes information about large 
deviations, which are encoded in the free energy function, defined for 
s E [-K, 1) - where k E (0, 1) is such that Av(V(0) K ) < oo - by 

M*;E) = lim **Mmof,E + iQ,»m (11 _ o) 

\x\— >oo \x\ 

and for s = 1 by (j) X (l; E) = lim^! 4> x (s; E). The existence of the limit in 
(11-a), for a.e. E E R, is shown in ( |AWb] . Theorem 3.2). We have 

Theorem 4c (Theorem 2.5 of [AWbj) Under the previous assumptions 
(but, much more generally, assumptions A-E, pg.8, of |AWb] ) . for any A > 
and Lebesgue a.e. E E R at which 

cj )x {l-E)<-\ogK (11-6) 

the Green's function satisfies almost surely 

9fG A (0, 0; E + iO, w) > (11 - c) 

By convexity arguments (see, again, theorem 3.2 of [AYVbJ), cj) X (s;E) > 
—sL x (E), and, hence, (11-b) implies (10-a) and, therefore, Theorem 4a is 
satisfied under the assumption of theorem 4a. 



(10 - d) 
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The importance of the free energy function is two-fold: it characterizes 
almost completely the different spectral regimes (see jAWbj ). and it is also a 
powerful tool in the dynamics, a topic to which we now turn our attention. 

2c - Dynamics 

We have seen the importance of the dynamics for spectral theory already 
in section 2a, through quantities like the sojourn time (5) and the idea of the 
proof of theorem 2. This is also true regarding theorem 3, which makes use 
of the fact that the FS transform of the spectral measure (associated to the 
Kronecker vector Sq for instance), defined as in proposition 1, satisfies 



by the Plancherel theorem and Fatou's lemma (see |Kle98] or exercise 3.5 of 
|MW12bj ). Klein succeeded in proving that 



with probability one for certain values of (E, A) mentioned in the statement 
of theorem 3, from which the spectral measure restricted to (— E, E) is purely 
a.c. by proposition 1. 

Another interesting dynamical consequence of (11-c) was emphasized by 
Miller and Derrida |MD94j . in a basic theoretical paper on weak disorder 
expansions: if quantum particles are sent coherently at energy E = k 2 + 
U w j re down a wire attached to the graph at the vertex x, the reflection 
coefficient Re satisfies Re < 1 iff (11-c) holds, see also the discussion in 



A different class of problems involving the dynamics relates to quan- 
tum dynamical stability. We briefly summarize the main issues, follow- 
ing ( |MW12b] . section 3d). Most of the concepts of dynamical stability in 
quantum mechanics are related to the growth of expectation values 



of certain observables A = A(0) as time t evolves. Boundedness of (A(t)) 
is related to quantum stability. In general, for a lattice model such as the the 



< — lim inf 

7T »?->-0+ 





[AWllbj . pg.3. 



(At)) 



(12 -a) 
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Anderson model if the potential is random, or the case of an almost-periodic 
potential, it is common to adopt A = \X\ 2 , where X denotes the operator of 
multiplication by the coordinate X, and one distinguishes the localization 
regime in which (|X| 2 ) < C for all t, or the ballistic regime where 

(|X| 2 ) = 0{t 2 ) , ast^oo (12-6) 

which are supposed to correspond, respectively, to p.p. and a.c. spectrum; 
Intermediate regimes are characterized by the diffusion constant 

= lim (12 - c) 

t— >oo t 

when the above limit exists , and the dependence on the initial state \I/ is 
explicitly indicated. 

If (|X| 2 ) < C with C a positive constant for all t, then \I/ has no continu- 
ous component. This is corollary 2.3.1. of [Las96j, and implies that = 0. 
It is remarkable that the converse, i.e., that if \I/ G "H pp , then = 0, 
is not true. Indeed, the only general result about the converse is Simon's 
paper |Sim90] on the absence of ballistic motion, i.e., lim t _ > . 00 (|.X'| 2 (£)) = 
if ^ G "Hpp, which is far from the expected = 0. The reason for this 
is the instability of the two exotic spectra, dense p.p. and s.c; indeed, 
even a rank one perturbation with arbitrarily small norm is able to induce 
a transition from one type to the other! ( |How92(, ISW86] ) . Since one does 
not expect the dynamics to be strongly affected by such perturbations, it 
is plausible that the absence of ballistic motion - the latter characterizing 
a.c. spectra - is both a feature of s.c. spectra and the (thick) p.p. spec- 
tra obtained from s.c. spectra by such perturbations: therefore Simon's 
result might be optimal! This has indeed be shown in the remarkable paper 
[dRJLS96j, where a potential was constructed such that the Hamiltonian H 
on Z 2 (Z) has a complete set of exponentially decaying eigenfunctions but, for 
any 5 > 0, | \X exp(— itH)8 \ \ 2 /t 2 ~ 5 is unbounded as t — » ±oo. Note that 
||xexp(— itH)5o\\ 2 = (exp(— itH)5o, X 2 exp(— itH)5o). In the words of Jito- 
mirskaya |Jit07] . this example showed that mere "exponential localization" 
of eigenfunctions need not have any consequence for the dynamics. Thus, 
[dRJLS96j was pioneer in demonstrating the importance of dynamical lo- 
calization , which, since then, was proved for various random models in the 
form 

suplHXl 9 Ei(H u ) exp(-itH u ty)\\ < oo 
t 
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with probability one. I C R is an energy interval in the localization region 
and Ej are spectral projections for /. For the Anderson Hamiltonian, or the 
deterministic almost-Mathieu model, several results are known, and we refer 
to ( |Jit07j . p. 621) for a discussion and (numerous) references. 

Concerning now the diffusion constant given by (12-b), one may define 
the so-called diffusion exponents (see, e.g., |BSB00j ): 



± log(T m ) 

where lim-t = limsup T _ 5 . 00 (liminfT- 5 . o)- Both (3 m {^) are nonincreasing func- 
tions of m and obey 

< /?-(*) < /?+(*) < 1 

(12 - d) 

The following important theorem ( |Gua89j . |Las96] . |Com93] ) relates the 
above quantities to the Hausdorff dimension of the s.c. spectrum: 

Guarneri-Last-Combes (GLC) theorem If H is self-adjoint on % 
with spectral family {E(\} and uniformly alpha Hoelder continuous spec- 
tral measure /jl^(X) = (^f, E(\)ty) (definition 1) there exists a constant c* im 
depending on \1/ and m such that, for all T > and m > 0, 

((*(t), \X\ m ^(t))) T > c^ m T ma / d (12 - e) 

Above, the time average < . >t is defined as in the statement of the 
Strichartz-Last theorem. By the GLC theorem, 

/?-(*) > a/d (12 - /) 

For fiis, a.c, a = 1; hence, in dimension d = 1, a.c. spectrum implies 
ballistic transport, while in dimension d = 3, a.c. spectrum and subdiffusive 
motion (/3 m < 1/2) may coexist. This motion is experimentally found in 
the so-called quasicrystals (see [GSB99J for a discussion and references). 
One of the very few models in d = 3 dimensions exhibiting a.c. spectrum 
with subdiffusive motion this was constructed in [BSB00J along lines similar 
to those of theorem 2, but using one- dimensional Jacobi matrices with self- 
similar s.c. spectra and making use of the aforementioned theorem of Hu 
and Taylor |HT94] . 
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Defining the transition probability 

P %t {x) = \(exp(-itH)^)(x)\ 2 (13 - a) 

a different form of the time- average (for the observable \x >< x\) is often 
easier to study |AWa] : 



POO 

P^^(x) = 2r] exp(—2r]t)Pq,j(x)dt 
Jo 

J | ((H - E- irjY 1 ^) (x) \ 2 dE 



1 

7T „ 

(13 - b) 

with inverse time parameter rj > playing the role of 1/T in the aforemen- 
tioned time- averages: the limit of long time averages T — > oo becomes 77 \ 0. 
The equality (13-b) follows from the spectral theorem and Plancherel's the- 
orem and explains the utility of this type of averaging. The averages (12-a) 
may now be written, with A(0) = \Xf, 



M^(3,t) = J2\^ P *, 



x 



and the modified time-averages as 

Ballistic (resp. diffusive) transport (12b)(resp. (12c)) correspond to r = 1 
(resp. r = 1/2) in the formulas: 

My(f3,t) ~fP (14 -a) 

and 

M 9 {P,V) ~V~ rP (14-6) 

in the finite-dimensional case Q = Z d . 
A general ballistic upper bound 

M^P,V)<C(P,^) V -^ (14 -c) 
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for all normalized ^ G 1 2 {Q) with XLee M^l^t^)! 2 < °°; follows from general 
arguments related to the Lieb- Robinson bounds |NS06j . see the appendix B 
of [Ma]- 

The strategy of proving the lower bound of the same type of the r.h.s. 
of (14-c) (with the purpose of showing ballistic behavior) is of particular 
interest. By the GLC theorem, for Q = Z d , in the case of uniformly alpha 
Hoelder continuous spectral measures, the Guarneri bound (see also |Gua 89j) 
becomes 

P*, v {x) < Cif (15 - a) 

, and thus 

Pr %v (\x\ < br ald ) = Kv&'t) 

\x\<br]- a / d 

<C d b d 

(15 - b) 

which implies that (14-a,b) can only hold if r > a/d, because, selecting b 
such that Cdb d = 1/2 one obtains My([3,r]) > ^ n ^ 1 . This lower bound 
diminishes, however, as d — > oo, and provides no information in this limit 
(which is the tree graph case). A different method becomes, therefore, nec- 
essary, of which we sketch the main ideas. With the notion of a.c. spectrum 
defined by 

a ac (H) = {E G R|Pr(3G(0, 0; E + iO) ^ 0) > 0} (15 - c) 

where Pr denotes the probability associated to the random potential (we 
refer to [AWaj for further discussion of this concept, but remark that a ac (H) 
is a nonrandom set which is a.s. the support of the a.c. component of the 
spectrum), the following theorem holds: 

Theorem 5 (Theorem 1.1 of |AWaj ) For any initial state of the form 
^ = f(H)5o, with a measurable function / G £ 2 (R) supported in a ac (H), 
and all b > 0, 

Av{Pr %ri (\x\ < br]- 1 ) < C(f)b + o{t]) (15 - d) 

with some C(f) < oo and o(t]) a quantity which vanishes for 77 — ^ 0. 

Theorem 5 yields a lower bound to M^,(/3, rj) of the ballistic form r.h.s. of 
(14-c); together with (14-c) it implies ballistic behavior for the full regime of 
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the a.c. spectrum. Its significance is due to the following facts: i.) in the hy- 
perbolic geometry of a regular tree, the classical diffusion (which corresponds 
to r = 1/2 in (14) in the finite dimensional case) also spreads ballistically, 
since, at each instant, there are more directions at which \x\ would increase 
than the one direction at which it goes down ; ii.) it excludes the possibil- 
ity of a different behavior in the previously discussed regime where the a.c. 
spectrum is caused by rare resonances, such as around the edges, where the 
density of states is very low and Lifschitz tail asymptotics hold (see |AWllaj ). 

The following diffusive bound is (for a bounded measurable I C a ac {H)) 
the main ingredient of the proof of theorem 5: 

C(I) 

esssup C( ; I+i{ o,i]Av{\G{0,x;()\ 2 ) < (16 - a) 

where C(I) < oo. The r.h.s. of (16a) coincides with the mean value of the 
total time spent at vertex x for a particle undergoing diffusion at the root 
x = 0. This bound is a consequence of the following bounds for the free 
energy function (11a): 

C) < -s log K for all ( G C + 
and for all s G [0,2] 

(16 - b) 

and 

Av(\G(0,x-()\ s < C + ( S ,C)exp[0( S ;C)|x|] 
for all s G [0, 1) and x G Q 

(16 -c) 

where C + (s, () G (0, oo) uniformly in ( G [-E, E]+i(Q, 1] for any < E < oo 
at a fixed s G [0, 1). 

Since for energies in the p.p. spectrum Av(\G(0, x; E + i0)\) = oo (see, 
e.g., pak06], theorem 3.15), (16c) cannot be expected to hold beyond s — 1, 
but theorem II. 2 of |AWaj shows that this bound extends to all energies such 
that 

esssup^+m^Avi^GiO, 0; 0)^) < oo (16 - d) 

for some 5 > 0, and then theorem II. 4 of |AWaj shows that (16-d) includes 
any bounded measurable set / C cr ac (H), with <r ac (H) given by (15-c). 
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2d - Concluding remarks on the Anderson transition 

An essential point, noted in |Jit07] . is that Klein's proof of theorem 3 
uses the looples character of the graph in an essential way and, thus, the 
Bethe lattice, while corresponding to the limit d — > oo, is, in a sense, one- 
dimensional. Similarly, the sparse model of section 2bl also has no loops, 
because theorem 2 "inherits" the one-dimensional structure of the model of 
theorem 1. Thus, both models still run short of solving the central issues of 
the original model, which remains one of the basic challenges of mathematical 
physics. 

The model on the Bethe lattice describes the limit d — > oo and, as re- 
marked in |AWllb] , the new and surprising phenomena found there may even 
have wider applicability due to the analogies drawn between tree graphs and 
many-particle configurations [BAL97J. They are, however, not expected in 
finite dimensions, in particular the absence of the mobility edge for weak 
disorder (Theorem 4b) is not found in the sparse model (see theorem 1 and 
theorem 2, where the existence of p.p. spectrum in a nonempty energy in- 
terval for arbitrarily small disorder is pictured as a concrete example of the 
tunneling instability). However, the methods used to prove theorem 4b, in 
particular concerning the free energy function (11a), may have wider applica- 
bility, in particular to investigate the existence or not of a sharp mobility edge 
in other models, such as the one in section 2bl or Krishna's model |Kricj . In 
addition, the ideas developed for the Bethe lattice dynamics reviewed in 2c, 
such as related to the negative moments of the Green's function (16-d), also 
deserve special attention as potentially powerful methods to prove existence 
of the a.c. spectrum. 

On the other hand, sparse models such as the one depicted in section 2bl, 
in dimension d > 2, comprise a part of the full model on a Bravais lattice. 
Although exponential sparsity is too severe, it should be recalled that the 
separable model (6) does not take account of dimensionality in a proper way, 
because in "truly d- dimensional" sparse models the cardinality of the set 
in (4) may change from O(log-R) to 0(R d ~ e ), for some e > 0, and "filling 
in" the remaining points may be expected to remedy this flaw, but this is 
an immensely challenging open problem! The sparse version has, however, 
one promising feature as compared with the full Anderson model: in the 
latter, the d > 2 version built as in (6) from the one- dimensional version (the 
separable full Anderson model) continues to have purely p.p. spectrum, by 
[G MP77| and the fact that convolutions of p.p. measures are p.p., in complete 
disagreement with the expected transition. It may thus seem worthwhile 
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to study sparse models such as that of 2b 1 as prospective descriptions of 
lightly doped semiconductors, which are also expected to exhibit an Anderson 
transition in dimension d > 2 [SE84]. Of course, non-sparse models such as 
the one proposed in |Kric] . are also very promising, although the available 
picture of the transition seems to be, as yet, somewhat less detailed than the 
model described in section 2bl. 

3-Frustration and short-range spin glasses: the Edwards- Anderson 
model 

Dilute solutions of atoms of large magnetic moments (such as the tran- 
sition metals Fe, Co, Mn) in a paramagnetic substrate (Cu, Au) present a 
number of peculiar physical properties. For small, but sufficiently high con- 
centrations of magnetic impurities, the susceptibility in low fields displays a 
characteristic peak, with discontinuous derivatives, at a temperature T g . The 
specific heat is always smooth, with a linear dependence on the temperature 
as T — > 0. The behaviour of these spin glasses has been explained in terms 
of an indirect RKKY interaction between the spins of these magnetic impu- 
rities mediated by the electrons of the paramegnetic matrix |B Y86j . which 
is of the form —J(\i — j\)SiSj with J(r) = (kpr)~ 3 cos(2kpr), where r is the 
distance between magnetic atoms and hp the fermi momentum; Si = ±1 
are (e.g.) Ising spins. The rapid oscillations and the weak decay of J(r), 
as well as the random distribution of the impurity magnetic atoms, are the 
basic ingredients of spin glasses. At sufficiently low temperatures there is 
a "freezing" of the magnetic moments in random directions (which leads to 
an increase of the susceptibility). The spin glass phase may be regarded 
as a conglomerate of blocks of spins, each block with its own characteristic 
orientation, in such a way that there is no macroscopic magnetic moment 
[BY86] . 

Edwards and Anderson |EA75j proposed a spin Hamiltonian, which we 
write in a generalized version as follows. For each finite set of points A C Z d , 
where the dimension d will be restricted to the values d = 2 and d = 3, 
consider the Hamiltonian 

Ha({J}) = £ 4 A* (17 -a) 

ijeA 

where <&ij with i,j G A are self-adjoint elements of the algebra generated by 
the set of spin operators, the Pauli matrices of, of , of, i G A, on the Hilbert 
space Ha = ®j e ACf , given by 

$ y = a x afa] + a y o\o) + a z a-a Z j (17 - b) 
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for = 1, and zero otherwise. The random couplings Jij, with = 1, 
are random variables (r.v.) on a probability space (Q,J-,P), where J 7 is a 
sigma algebra on Q and P is a probability measure on (O, J 7 ). We may take 
without loss of generality 

tt=x Bd S (18) 

where S is a Borel subset of R, B d is the set of bonds in d dimensions, and 
assume that the Jij are independent, identically distributed r.v.. In this 
case, P is the product measure 

dP=x Bd dP (19) 

of the common distribution Pq of the random variables, which will be denoted 
collectively by J. The corresponding expectation (integral with respect to 
P) will be denoted by the symbol Av. We have to assume that 

Av(J itj ) = (20) 

for all i,j G Z d , i.e., that the couplings are centered. This assumption 
mimicks the rapid oscillations of the RKKY interaction. Let E\ denote the 
GS energy of i.e., E\ = inf spec(H\). The following result was proved, 
among several others, in |CL10j : 

Theorem 6 ( |CL10j ) For P- a.e. {J}, the limit below exists and is 
independent of the b.c: 

e W = l im ^ (21 - a) 

Ayz<i |A| ' 

and 

e (d) = i n f (21 - b) 

a |A| v ' 

where |A| denotes the number of sites in A. Finally, 

e (d) > e (d+i) _ c ) 

(21-a) is the far-reaching property of self-averaging (see | And 78] for a dis- 
cussion): it expresses that measurable - e.g. thermodynamic- quantities are 
the same for any typical configuration of the sample, i.e., are experimentally 
reproducible. It follows from (21a) that, P- a.e., 

C M = lim ^1 (21 - d ) 

A/*Z<i |A| v ; 
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Let An denote a square with N sites if d = 2 or a cube with N sites if d = 3, 
and write if# } (.7) = H Aff ( J )- We 

now adopt periodic b.c. for simplicity, but 
the final result is independent of the b.c. due to theorem 1. We may write 

H${J)= £ H<f>{J) (22 -a) 

where Hn^ is given by 

H^(J) = c d J i,&J ( 22 

Above, q are factors which eliminate the multiple counting of bonds, i.e., 

c 2 = 1/2 (22 - c) 

and 

c 3 = 1/4 (22 - d) 

(2) 

and An is a square labelled by a site n, for which we adopt the convention, 
using a right-handed (x,y) coordinate system, that n = (n x ,n v ) is the vertex 

(3) 

in the square with the smallest values of n x and n y . Similarly, A„ is a 
cube labelled by a site n = (n x ,n y ,n z ) with, by the same convention, the 
smallest values of n x , n y and n z . Due to the periodic b.c, the sum in (22- 
b) contains precisely N lattice sites. The notation is short-hand for 
its tensor product with the identity at the complementary lattice sites in 
A N \A { n d) . Let E { N ' d \j) denote the GS energy of H {d) (J) , E^' d) (J) the GS 
energy of Hn\j) and 

E [d) =Av(E^ d) {J)) (22 -e) 

By the condition of identical distribution of the r.v. Jjj, E^ 1 does not depend 
on n, which is implicit in the notation used. We have 

Theorem 7 (Theorem 2 and proposition 1 of |Wrel2j ) The following 
lower bound holds: 

ME<r\j)) £ E ^ (23 _ a) 

Further, let dPo in (9) be the Bernoulli distribution dP = 1/2(5 j + 5-j), set 
for simplicity J = 1 and let a x = a y = in (17b). Then, 

,(2) 



e {2) > -3/2 (23 - b) 
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and 

e <3> > _ 1/4 36096 

' 4096 v 1 

The special case a x = a y = in (17-b) is the classical EA spin-glass; we 
also set a z = 1. In this case given by the r.h.s. of (21-d), is invariant 
under the "local gauge transformation" erf — > — of together with J i; j — > 
—Ji,j for all j\\j — i\ — 1, whatever the lattice site i. 

According to the above, |Tou77] . an elementary square ("plaquette") P 
is said to be frustrated (resp. non-frustrated) if Gp = Yl^^eP^d = ~^ 
resp. +1. Note that Gp is gauge- invariant and that, for the quantum XY 
(or XZ) model defined by setting a y = in (17-b), is also locally gauge- 
invariant if we add to the above definition the transformation erf — > —erf, 
i.e., the transformation on the spins is defined to be a rotation of 7r around 
the y - axis in spin space. The property of local gauge-invariance guarantees 
the absence of a macroscopic magnetic moment (spontaneous magnetization) 
mentioned before as a basic property of spin glasses |JAS81j . 

Since the Pauli z-matrices commute, finding the minimum eigenvalue of 
(17-a) in the classical EA case is equivalent to find the configuration of Ising 
spins c>i = ±1, denoted collectively by a, which minimizes the functional 

F(er, J) = JijCTiCTj 

a 

The minimal energy of a frustrated plaquette equals Epj = —2 and of a 
non-frustrated plaquette Ep tU = —4. 

(23-b,c) of theorem 7 provide the first (nontrivial) rigorous lower bound 
both for and e^ 3 -*. Using the natural misfit parameter 

m = J — ,' ' 01 24 
\E ld \ K 1 

as a measure of plaquettes frustrated or bonds unsatisfied (see (4) of |KK95j ). 
where E denotes the ground state energy of the frustrated system and E ld 
is the ground state energy of a relevant unfrustrated reference system, we 
find from (23-b) in the d = 2 case the lower bound m > 0.25 and for d = 3, 
from (23-c), the lower bound m > 0.26 • • • : thus, in both cases, the measure 
of frustrated plaquettes or unsatisfied bonds as defined above is at least of 
the order of 1/4. 
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The method of proof of theorem 7 is a rigorous version of a finite-size 
cluster method, originally due to [BS79J, together with the variational prin- 
ciple. If we take in (17) a z — 1, a y — 0, and a x = a, and consider a as a 
small parameter, we have the anisotropic XZ (or XY) model. By the norm- 
equicontinuity (in the volume) of H\({J}) (given by (17)) as a function of 
a, which is preserved upon taking averages over the probability distribution 
of the J, it follows that (23-b,c) hold with the right hand sides varying by 
small amounts if \a\ is sufficiently small: this is conceptually important for 
reasons mentioned in the introduction. 

The mean field theory, recently rigorously solved (see |Gue06j for a re- 
view) does not exhibit frustration - indeed, this concept is not even generally 
defined for their model, since the theory does not require an underlying lat- 
tice. Whether frustration is an important issue in the description of realistic 
spin glasses is an important open problem. 

We refer to the article by Bovier and Frohlich [BF86J (see also Bovier's 
book |Bov06j ) for an illuminating discussion of complementary, mostly global 
(i.e., involving the lattice as a whole) aspects of frustration. Our bound 
(23-c), which relies only on the local structure, is, however, slightly better 
than Kirkpatrick's [Ki r77j . which is based on very reasonable, but unproved 
conjectures of a global nature. A most relevant open problem would be, of 
course, to extend the finite-size cluster method to obtain bounds for the free 
energy of the EA model. 

4- Return to equilibrium and (thermodynamic) instability of the 
first and second kinds 

The dynamics of spin glasses is a topic of great relevance, both conceptual 
and experimental. It happens, however, that the standard approach to the 
kinetic theory (e.g., for the mean-field spin glass model) relies on Glauber 
dynamics (see, e.g., |Sza97] ).a well-known dynamics imposed on the Ising 
model, which has been used to study metastability in early days |DCQ74] and 
more recently |SS98] . In spite of the considerable independent mathematical 
interest of the latter works, it should be remarked that only for a quantum 
system does a physically satisfactory definition of the dynamics of the states 
and observables exist which is relevant to the microscopic domain, including, 
of course, condensed matter physics. The scarcity of examples of approach to 
equilibrium in quantum mechanics is due, of course, to the extreme difficulty 
of estimating specific properties of the quantum evolution. 

This fact was the basic motivation which prompted Emch |Emc66j and 
Radin |Rad70] to propose a quantum dynamical model (which we dubbed 
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the Emch-Radin model in |Wrel2j and |MW1 2b|). which is relevant to sys- 
tems with high anisotropy, and displays a remarkable property of non- 
Markovian approach to equilibrium, or return to equilibrium. It turns 
out that such a model is useful not only for a description of nuclear spin- 
resonance experiments - such as the one |LN57j which motivated Emch and is 
described below - but also to describe dynamical effects associated to quan- 
tum crossover phenomena in spin glasses (see |Sac94] for a review and ref- 
erences). For this purpose, a material was chosen with a strong spin-orbit 
coupling between the spins (of the magnetic ion) and the underlying crystal: 
this coupling essentially restricts the spins to orient either parallel or antipar- 
allel to a specific crystalline axis, which we shall label as the z-axis. Such 
spins are usually referred to as Ising spins, and will be described by the Ising 
part of the forthcoming Hamiltonian. In the experiments (see [Sac94j and 
references), a transverse magnetic field is then applied, oriented perpendicu- 
lar to the z axis. A large enough transverse field will eventually destroy the 
spin glass order even at T = 0, leding to the existence of a crossover region. 

We shall be interested in a situation in which a large transverse field 
has been applied; the initial state of the system is, then, approximately, a 
product state of the forthcoming form (16). The same experimental setup 
should therefore allow a measurement of the rate of return to equilibrium 
of the mean transverse magnetization, i.e., of how f(t), defined by (20), tends 
to zero, or of how (21) is approached. As we shall see, this rate depends 
strongly on the probability distribution of the J. 

The XY model has also been extensively studied from the point of view 
of return to equilibrium, since early days as an example of general theory 
|Rob73] . more recently in [AB73j (see the latter for reference to related im- 
portant work of Araki), but unfortunately only in one dimension. 

The nonrandom model 

Following Radin's review |Rad70] of the work of Emch |Emc66] . consider 
the experiment of |LN57j : a CaFi crystal is placed in a magnetic field thus 
determinig the z-direction, and allowed to reach thermal equilibrium. A 
rf pulse is then applied which turns the net nuclear magnetization to the 
x direction. The magnetization in the x direction is then measured as a 
function of time. 

As in ( |Emc6 6j . [Rad70] ) . we assume an interaction of the form 




(25) 
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where j G V and k G V in (15), with V is a finite region, V G Z d . i7y 
is defined on the Hilbert space % = <S> iGV C 2 . The state representing the 
system after the application of the rf pulse will be assumed to be the product 
state 

P = P(O) = (g)0i (26 -a) 

jev 

where 



^•(•) = tr 3 (- exp(- 7 aj))/tr,(exp(- 7 aj) (26 - b) 

and tr is the trace on C 2 . Other choices of the state are possible |Rad70] . 
but we shall adopt (16) for definiteness. Let 

S' = l/N(y)Y,°j (27) 
jev 

be the mean transverse magnetization, with N(V) denoting the number of 
sites in V. The real number 7 in may be chosen as in [Emc66j such as to 
maximize the microscopic entropy subject to the constraint trS x p(0) equal 
to a constant, i.e., a given value of the mean transverse magnetization. Since 
the state (26) is a product state, 7 is independent of V, and has the same 
value if S x is replaced by af Q , for any i G V. Let 



and 



p v t = IfypiftU? (28 - a) 



<)v(t)= P n<) (28-6) 



where Uv(t) = exp(itHy). We may write, by (18), 

«)v(t) = p(0)(U v (-t)alU v (t)) (28 - c) 

It is natural to define 

«}(*)= lim«)v(t) (29) 

provided the limit on the r.h.s. of (29) exists, as the expectation value of 
the local transverse spin in the (time-dependent) nonequilibrium state 
of the infinite system. A weak-star limit of the sequence of states Py{-) = 
tr(-Uyp(Q)Uy ) exists, by compactness, on the usual quasilocal algebra A 
of observables (see, e.g. |Sew86j . for a comprehensive introduction to the 
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algebraic concepts employed here). The expectation value of erf in the equi- 
librium state associated to (25) is zero by the symmetry of rotation by 7r 
around the z axis. We may now pose the question whether the limit 

f(t)= hm pUVN(A)J2<) (3°) 

where A denotes a finite subset of Z", which is interpreted as the mean 
transverse magnetization, exists. The property of return to equilibrium is 
expressed by 

hm /(*) = (31) 

Of particular interest is the rate of return to equilibrium. In the present 
nonrandom case, the limit at the r.h.s. of (31) equals f(t) = p t oo (o'f ), for 
any i , by translation invariance of p^. This is not so for random systems, 
in which case additional arguments are necessary to show the convergence of 
the r.h.s. of (31) for almost all configurations of couplings. 
We assume that the function e(-) satisfies: 



and either: 



or 



e(0) = (32) 

^2 < oo i.e. e E h(Z d ) (33 - a) 

jez d 



5>(j)) 2 < oo i.e. eel 2 (Z d ) (33-6) 
jez d 

The random model 

We now introduce the Hamiltonian of a disordered system corresponding 
to (25), where we take, for simplicity, 

V = V n =[-n,n} d (33-d) 



H n = 1/2 J UMJ - k )^t (34) 

where J<j,k) are independent, identically distributed random variables (i.i.d. 
r.v.) on a probability space which we denote by (fl,B,P). We shall use 
Av(-) to denote averaging with respect to the random configuration {Jj,k}, 
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which we denote collectively by the symbol J. The J are assumed to satisfy 
|vEvH83j : 

Av(J m )=0 (35 -a) 

\Av{J$ >k) )\ <n\c n Vn = 2,3,4,- ■■ (35-6) 
We shall take, without loss, e(j) > 0. If 



<3) 



p if j e ±<5j , 

otherwise , 



(35 - c) 



where ±5j for z = 1, • • • , d denote the set of z = 2d bonds connecting the 
origin to a point of Z d , we have the EA spin glass model of section 3. 

Let \V n \ denote the number of sites in V n and the free energy per site /„ 
be defined by 

n / t\ -kT\ogZ n (J) 

fn{J) = prr, (36 - a) 

I n I 

where 

Z n (J) = tr(exp(-pH n )) (36 - b) 

is the partition function and the trace is over the Hilbert space "H. Then 
Theorem 8[vEvH83] (see also [KS79J for the first result in this direction): 
Under assumptions (22-b) and (35), the thermodynamic limit of the free 

energy per site 

f(J) = lim f n (J) (37 -a) 

n— >oo 

exists and equals its average: 

f(J) = Av(f(J)) = lim Av(f n (J)) (37 - b) 

n—>ao 

for almost all configurations J ( a.e. J). 

The reason why (33-b) suffices for the existence of the thermodynamic 
limit is that, in order to obtain a uniform lower bound for the average free 
energy per site, the cumulant expansion (see, e.g., ( |Sim79j . (12.14), pg 129, 
for the definition of Av c , there called Ursell functions): 



Av(exp(tJ (itj) ) = exp(^^ c (J^. ) )r/n!) 



n=2 
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was used |vEvH 83j . which, by (35-a), starts with the second cumulant 
Av c (J^ j\) = Av(J^^), which is the variance of J(ij). Condition (35-b) was 
used to control the sum in the exponent above. 
Thermodynamic stability 

Stability considerations play a key role in quantum mechanics [LS10J. 
Let a statistical mechanical system be described by a collection of amil- 
tonians associated to finite regions A C R d , for particle systems, or 
A C Z d for spin systems, with volume V^(A): examples are (25) and, for 
random systems, (34): in the latter case there is implicit in the Hamiltonians 
a dependence on the random variables J, and the constant c in definition 2 
below is assumed to be a.e. independent of J. The system's free energy is 
/ A = ~ fcT y( g A ^ A( ' J ' > with Z\ = tr(exp(— /3 H/Cj), and the thermodynamic limit 
means A/R d (or Z d ) with the proviso of fixed density for particle systems, 
where A /* Z d denotes a limit in the sense of van Hove or Fisher (the latter 
being required for random systems, see the discussion in |vEvH 83j : (33-d) 
satisfies this assumption). 

Definition 2 The system is said to be thermodynamically stable if 
there exists a constant < c < oo such that 

/a > -c (38 - a) 

It is said to satisfy thermodynamic stability of the first kind if 

H A > -oo (38 - b) 

and to satisfy thermodynamic stability of the second kind if 

H\ > -cV(A) (38 - c) 

The above definitions (38-b,c) are patterned after the corresponding ones 
for N- body systems in |LS10j : they also apply to relativistic quantum field 
theory, where the particle number N is not conserved. 

Theorem 5 illustrates the interesting fact that, for disordered systems, 
stability of the second kind (38-c) may fail even for a thermodynamically 
stable system: this happens when only (33-b) (but not (33-a)) holds. In 
the following, we shall demonstrate that this phenomenon has important 
dynamical implications for a relevant class of disordered systems. 

We now return to our nonrandom model (25). When (33-a) holds, it fol- 
lows from a folklore theorem (see, e.g., Theorem 6-1 of |MW12b] or Theorem 
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3.3, pg. Ill, of |Dav76j ) and a representation of / ( |Rad70| . pg.295, and 
proposition 1) that exponential decay in the sense that 



|/(f)l<Cexp(-d|f|) (39) 

with C and d positive constants, cannot hold. This is essentially the condi- 
tion that the physical (GNS) Hamiltonian Hq is bounded from below (semi- 
bounded) (see again |Sew86j ). where 

H = al Y, <kK (40 - a) 

The infinite sum (40) stands for a limit in norm in the quasi-local algebra A 
associated to the spin system, see [Rad70], loc. cit. It is the thermodynamic 
stability of the second kind (38-c) which guarantees the existence of Hq, and 
(39) as a consequence. In the random case (34), 

Ho(J) = < Yl J ^e{k)a% (40-6) 

under assumption (33-a): if only (33-b) holds, the representation (40-b) is, of 
course, not defined. We refer to jWre!2| . proposition 4.1 for the dynamical 
consequences of this fact: in particular, unlike the stable case (33-a), expo- 
nential decay (39) does indeed hold for a class of one-dimensional potentials 
with algebraic decay, as shown there. 

We now consider the special nearest-neighbor interaction (35-c) - i.e., the 
EA model of section 3. The following theorem may be proved (along the same 
lines, but simpler, than theorem 3.2 of |Wrel2j ) (There are two misprints in 
the latter reference: in the r.h.s. of (40a) the number four should be inside 
the exp function, while in (40b) the first parenthesis should be moved to the 
left of the symbol Av, and the last parenthesis omitted): 

Theorem 9 Let (35-c) hold, i.e., (34) describe the EA spin glass model 
of section 3. Then there exists a subsequence {m„} ne z such that, a.e. with 
respect to J, 

/(*) = J™ ..M^^-) 

* n 

SAv{ Po {t)) 

(41) 
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where 

5 = O «) ^ (42) 

and 

p (t) =Y[cos(2/3J 0>±Si t) (43) 

±Si 

We have the 

Corollary 9 Consider the distribution functions 

dPi(x) = l/2[5\(x) + 5-i(x)} Bernoulli distribution (44 — a) 

dP 2 /dx = l/2x[-i,i](x) uniform distribution (44 — b) 

dP^/dx = — jz — exp(— x 2 ) Gaussian of unit variance (44 — c) 

The corresponding values of /, defined by the r.h.s. of (41), are: 

f x (t) = (cos(2/3t) 2 (45 -a) 

/.(t) = (^)" (45-6) 
f 3 (t) = exp(-2zt 2 ) (45 -c) 

Remark 2 

a. ) In the Bernoulli case, (45-a) implies no decay, as in the non-random 
model: this is an instance of a result of Radin for general interactions sat- 
isfying (33-a), according to which / is almost-periodic ( [Rad70j . Lemma 2, 
pg.1951). 

b. ) For the uniform distribution, (35-b) describes a non-Markovian 
return to equilibrium. Clearly, in this case, under assumption (33-a), 
(40b) is a physical Hamiltonian which is semibounded for each realization 
of J, and thus exponential decay (39) cannot hold by the aforementioned 
folklore theorem. 

c. ) In the Gaussian case, not even thermodynamic stability of the first 
kind (38-b) holds, because the Gaussian r.v. range over R. (40b) shows that 
also the physical Hamiltonian Hq is not bounded below, and thus (39) may, 
in principle, hold. That this is actually the case is, of course, made clear by 
(45-c). 

d. ) Since the transverse magnetization in (41) is an observable (mea- 
surable quantity), self-averaging is an important requirement, as discussed 
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before. We conjecture that all subsequential limits are equal, a.e. J, to the 
r.h.s. of (41), i.e., that the limit n — > oo in (41) exists. 
The issue of probability distributions 

As is well-known, the use of Gaussian probability distributions (p.d.) in 
disordered systems is standard: it simplifies several passages considerably, 
e.g., the integration by parts formula in the first proof of the existence of the 
thermodynamic limit for the mean field theory in [GT02] . One might think 
that the results are not expected to depend qualitatively on the probability 
distribution: this is, however, not so. 

In |SW85j , where a mean- field Ising model in a random external field was 
studied, it was proved that the existence or not of a tricritical point in the 
model's phase diagram is tied to the probability distribution: it is absent for 
a Gaussian p.d., and present for a Bernoulli distribution. 

On the other hand, the results on the Anderson model on the Bethe 
lattice reviewed in 2b2, 2b3 pose absolute continuity of the probability dis- 
tribution as a major requirement. Also, as remarked by Jitomirskaya ( |Jit07j . 
pg. 624), "proving localization for the Anderson model with a Bernoulli dis- 
tribution remains more than a mere technical challenge" (this is problem 3 
in Jitomirskaya list of open problems). 

As an attempt to explain the result in [SW85j physically, we conjectured 
that it was due to the fact that a discrete distribution of probabilities such as 
the Bernoulli distribution samples just a few values of the couplings and thus 
introduces some short-ranged elements into the problem: from this point of 
view, discrete distributions may have a closer connection with real materials. 
This feature is shared by the uniform distribution, see remark 2 b.), but 
even here there is a remarkable dynamical difference between Bernoulli and 
uniform distributions with regard to decay, exhibited by (45-a) and (35-b). 
The latter yields, somewhat surprisingly, the same result in one dimension 
(z = 2) as obtained by Emch [Emc66j with an interaction of infinite range 
(e(|n|) = 2"'™'" 1 ) in the nonrandom model (or, alternatively, the random 
model with Bernoulli distribution and the same potential). Since the latter 
(algebraic slow decay with oscillations) seemed to lead to a good qualitative 
description of the nuclear spin resonance experiment |LN57] . the uniform 
distribution decay (45-b) might well be in qualitative agreement with the 
previously proposed spin glass experiment! 

In contrast to a discrete or uniform distribution, a Gaussian distribution 
samples many values of the couplings and works to reinforce the long-range 
nature of the interactions. For effectively short-range inteactions, this might 
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not be suitable. In any case, the Gaussian probability distribution leads to a 
spectacularly fast rate of return to equilibrium (45-c), and such sharp differ- 
ences, such as between (45-c) and (45-b), should be amenable to experiment. 
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